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Abstract 



The circadian clocks keeping time of day in many living organisms rely on self-sustained biochemical 
oscillations which can be entrained by external cues, such as light, to the 24-hour cycle induced by 
Earth rotation. However, environmental cues are unreliable due to the variability of habitats, weather 
conditions or cue-sensing mechanisms among individuals. A tempting hypothesis is that circadian 
clocks have evolved so as to be robust to fluctuations in daylight or other cues when entrained by 
the day/night cycle. 

To test this hypothesis, we analyze the synchronization behavior of weakly and periodically forced 
oscillators in terms of their phase response curve (PRC), which measures phase changes induced by 
a perturbation applied at different phases. Wc establish a general relationship between, on the one 
side, the robustness of key entrainment properties such as stability and phase shift and, on the 
other side, the shape of the PRC as characterized by a specific curvature or the existence of a dead 
zone. This result can be applied to computational models of circadian clocks where it accounts 
for the disparate robustness properties of various forcing schemes. Finally, the analysis of PRCs 
measured experimentally in several organisms strongly suggests a case of convergent evolution toward 
an optimal strategy for maintaining a clock that is accurate and robust to environmental fluctuations. 
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INTRODUCTION 

Circadian entrainment is the process by which a biological clock with a period of approximately 24 
hours is synchronized to environmental cycles associated with Earth rotation. A stable and precise 
phase relationship between internal and external times is vital for organisms that need to timely 
anticipate dawn or dusk in order to coordinate their physiology to diurnal environmental changes. 
For instance, a precise clock in phototrophic organisms like cyanobacteria or plants has been shown 
to optimize cell growth and fitness (jH-IsI). 

However, circadian clock precision is challenged by many sources of intrinsic and extrinsic vari- 
ability, which could dramatically affect key properties of endogenous biochemical oscillations. A 
growing interest has therefore been devoted to investigate how the period and amplitude of these 
oscillations are impacted by genetic mutations (0, lEl), molecular noise (@-[9l), or contextual variabihty 
(@, [13)) which led to advocate the existence of design principles ensuring robust oscillatory behavior 



in circadian clocks (IllHlSl). 

Nevertheless, a robust endogenous clock does not guarantee by itself a precise phase relationship 
with the day-night cycle, since the environmental cues associated with the diurnal changes are also 
highly fluctuating. The daylight intensity and quality sensed by an organism depend on various 
environmental factors such as meteorological conditions, shade habitats or, for marine organisms. 



the distance to sea surface and water turbidity (jl4j . Il5l ). In addition, variations of the behavior 



and light-sensing abilities of individuals can also alter the light signal reaching their core molecular 
clock. Thus, individuals of the same species living in different places or times can experience a wide 
spectrum of light intensities over several orders of magnitude, raising the question of whether and 
how specific robustness strategies are implemented in their clock architectures to maintain a precise 
synchronization despite unreliable environmental cues. Although the importance of this problem was 
noted some time ago (flil), the robustness of circadian clocks to daylight fluctuations and how this 



constraint shapes their molecular architecture have been little studied until very recently (1171. Il8l). 
These computational studies have revealed disparate robustness capabilities depending on the clock 
model, which remains to be explained in a comprehensive approach. 

A natural theoretical framework to address this issue is the dynamical system theory of synchro- 



nization (Il9l. |20|). In forced oscillatory systems, a phase-locked state specified by a stable phase shift 
between the external and the internal phases closely depends on forcing features such as amplitude, 
profile or period. Changing these forcing properties can alter the phase shift or induce complex 
dynamical regimes (e.g., quasi-periodic, period-doubled, or chaotic), which are both undesirable for 
a functional circadian clock. Thus, the question of whether a circadian clock is robust to daylight 
fluctuations can be generalized to the problem of how an oscillator can maintain a stable synchro- 
nization and phase relationship with a forcing cycle that exhibits significant variability. Assuming 
a weak forcing allows one to drastically simplify the theoretical analysis by utilizing perturbation 



theory in the vicinity of a periodic orbit (I2ll423[). in which the infinitesimal phase response curve 



(IPRC) of a weakly coupled oscillator provides important information on its synchronization. 

In this paper, we first lay out the theoretical approach within which two quantities measuring 
robustness of the entrained state with respect to forcing fluctuations are defined. The phase approx- 
imation in the weak forcing limit allows us to identify the geometric properties of phase response 
curves that contribute to the robustness and precision of the clock phase in presence of fluctuations 
in the forcing amplitude. The general criteria obtained are shown to explain and predict the robust- 
ness properties of biologically-based circadian models. Finally, the analysis of PRCs that have been 
measured experimentally in several organisms supports the idea that living organisms have evolved 
salient strategies to maintain an accurate clock that is robust to daylight fluctuations. 
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RESULTS 

Circadian clock entrainment and phase approximation 

The entrainment of circadian clocks by cyclic environmental changes is a paradigmatic example of 
an unidirectional sjTichronization process. A biochemical implementation of this process in a living 
cell involves a network of genes and proteins interacting with each other, whose temporal evolution 
is usually modeled by a set of ordinary differential equations: 

^ = F(X(t),po + dp(t)) (1) 

where the components of vector X are the concentrations of the molecular actors interacting according 
to the biochemical kinetics F, derived from the law of mass action. 

Synchronization of the circadian clock to the diurnal cycle requires that some parameters pi are 
sensitive to temporal changes of light induced by the cycle, so that they differ from their value in 
the dark (po)i by 

dpi{t) = eL{u)vi{po)i (2) 

where the perceived light intensity is described by an amplitude e and an a normalized temporal 
profile L which depends on u = t — id, with u G [0,T], id and T corresponding to dawn time and 
the 24-hour day length, respectively. The relative sensitivities of parameters to light are given by 
the vector v, whose is normalized (||v|| = 1) to ensure that e is a relative modulation amplitude. In 
day/night entrainment conditions, the light sensed by the organism is assumed to be restricted to 
daytime of duration td {L{u) = for m G [td^T]). In this work, an important point is that e and L 
can differ between individuals depending on environmental or physiological context. 

The existence of endogenous circadian oscillations requires that Eqs. [1] parameterized by po in the 
absence of light have an asymptotically stable limit cycle solution Xq characterized by a free-running 
period (FRP) Tq. The light-dependent perturbation deviates the circadian oscillator from its free 
limit cycle trajectory during daytime. If the amplitude deviation is not too large, the stability of the 
limit cycle ensures that the deviation decays in the absence of perturbation (during night), with the 
only memory of the past perturbation being a residual phase change. Thus, the Poincare map of the 
dynamical system described by Eqs. [1] can be approximated by the following unidimensional map 



0„+l = G{(t)n) = 0n - 7 + ^(0n, e) (3) 



where (j)n G [0, T] is the oscillator phase in units of circadian hours (abbreviated as ch) (|24| ) at which 
the light-dependent perturbation is switched on during the n-th day (i.e., the oscillator phase at 
dawn). 7 = (Tq — T)T/Tq is the phase difference associated with the period mismatch between the 
forcing and endogenous periods. V{(f)n, e) is the phase change induced by the perturbation given by 
Eg. [2] applied at oscillator phase 0„, and is known as a phase response curve (PRC) in the litterature 



(1201 ). The mapping G has a stable fixed point if there exists a phase shift 0* that satisfies: 

-2<9<^r(0*,e) <0 

This phase-locked state is usually termed 1 : 1 synchronization (or entrainment) state. The stability 
coefficient d^V{(f)*,e) (also noted V'{(f)*)) characterizes how fast fluctuations around 0* decay. 

Due to the transversal stability of the limit cycle and the neutral stability of the phase, one can 
assume that small perturbations mainly induce phase deviations, which comprehensively describe the 
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effect on oscillator dynamics. In the framework of the so-called phase approximation (22), the PRC 
V{(f)) can then be derived from an impulse infinitesimal PRC (i-IPRC) Z{u) according to a convolution 



integral, with the important property that it scales linearly with light stimulus amplitude (jlll ): 



V{(f),e) = e Z{u + (j))L{u)du = eW{(f)) (5) 
Jo 

Besides the classic i-IPRC Z{(j)), Eq. [5] also introduces the d-IPRC W{(j)) that describes the phase 
response of the oscillator to day-like light perturbations of duration td, temporal profile L{u) and 
vanishingly small amplitude e. 

Two distinct metrics to measure clock entrainment robustness 

The daylight intensities sensed by a living organism depend on environmental factors, habitat or 
cloud cover, which can vary in time and space. They can also vary due to the temporal or inter- 
individual variability of the light-sensing mechanism. To determine the impact of this variability 
on the clock precision, our approach is to investigate how the PRC V{(f)) given by Eq. [3] reacts to 
changes of the daylight driving cycle, as recapitulated in Fig. [H Indeed, changes in amplitude e 
and temporal profile L{u) modify the asymptotic entrainment state (Fig. [H^), and the resulting 
variations in phase shift or stability coefficient can be determined from the associated changes in the 
PRC (Fig. mB). 

To define relevant measures of robustness, we consider small changes of daylight intensities with 
respect to some average value: 

eLiu) = eo{Lo{u) + 5eL{u)) (6) 

where u E [0, To], L is an appropriately normalized perturbation profile and 6e is small compared to 1. 
Such a perturbation generally modifies the PRC V{(f), e), and consequently the phase shift 0* and the 
stability coefficient V {</)*). Considering that these two important properties of entrainment should 
vary as little as possible, two complementary sensitivity measures 11 and S can be defined, which 
are respectively the variance of the phase shift and the relative variance of the stability coefficient 
in response to small daylight fluctuations of randomly distributed amplitudes (see Appendix ^4). 
Assuming that the phase shift and stability coefficient depend linearly on amplitude perturbation 6e, 
n and S can be written: 



(7) 



n = [D5e0*(eo)]' 
S = [Ds,d^V{(j)\eo)/d^V{(j)\eo)] 

The quantity S provides indirect information about the entrainment range of the forcing mech- 
anism. Indeed, there are usually minimum and maximum values of the forcing strength such that 
outside the interval they define, complex dynamical regimes such as quasiperiodic, period-doubled 



or chaotic ones occur (Il9l ). These alternative regimes are associated with a high variability of the 
phase shift, which is incompatible with the time-keeping function of the clock. Achieving stable 
entrainment for a large spectrum of daylight intensities therefore requires that the stability coeffi- 
cient is sufficiently independent of forcing strength so that it does not approach marginal stability 
{V'{(j)*) = or —2) too closely. Even when the 1:1 entrainment state is stable over a wide range, the 
phase shift 0* can still vary in presence of fluctuations of daylight amplitude and profile, which is 
measured by 11. It is important to emphasize that, because 11 and E are defined in with respect to the 
asymptotic phase shift, this approach is restricted to the case where the relaxation time toward the 
asymptotic state is smaller that the correlation time of daylight fluctuations over successive days, for 
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instance when individuals are subjected to slow changes of weather, habitat, genotype or phenotype. 
This approach may nevertheless give reliable insights into the case where daylight intensities change 
from day to day. 

n and E sensitivities depend on the shape of IPRCs 

The sensitivity measures defined by Eq. [7] characterize variations in phase shift and stability coef- 
ficients. These variations depend themselves on the response of the PRC V{(j)) to small changes in 
the amplitude and temporal profile of the forcing scheme, which must be computed numerically or 
determined experimentally. In the limit of weak forcing, however, the PRC can be approximated 
as a mere convolution of Z with the forcing profile (see Eq. [5]). Subsequently, S and 11 can be 
expressed in a simple manner as a function of forcing amplitude e, period mismatch 7 and IPRC 
characteristics. In the following, we investigate the robustness with respect to two types of dayhght 
fluctuations which may occur in nature, namely changes (i) in the daylight temporal profile or (ii) 
in the average dayhght amplitude. 

In the first case, the primary effect of the temporal profile changes in light signal is to produce 
phase-shift variability. 11 can then be expressed as the function of Z{(f)) (see Appendix B): 

_ r £° Z{u + (j)*)L{u)du 
~ [j^°Z'{u + ^*)Lo{u) du 

Low values of 11 are achieved by minimizing the ratio between the two integrals, which constraints 
function Z. The numerator of 11 in Eq. |8]can be minimized for arbitrary light input profile fluctuations 
L{u) when Z is constant over the coupling interval [(f)*, 0* -f- td], since J L{u)du = 0. In contrast to 
this, maximizing the denominator of 11 in Eq. E] typically requires that Z decreases significantly over 
+ T£,]. An optimal compromise between these two requirements is reached when the i-IPRC 
Z is constant on a large portion of the coupling interval except at the beginning and the end of 
the interval, where phase advances or delays occur, so that the dead zone slightly smaller than the 
coupling interval. To illustrate this point. Figs. [2JA. and [2p show two i-IPRCs computed for two 
different coupling schemes of the circadian model introduced in the following section. The IPRC 
in Fig. [2] A displays a large dead zone during daytime contrary to the one in Fig. [2] C. We will see 
later that this is associated with very different degrees of robustness of the two coupling mechanisms. 
Besides 11, quantities 11^ characterizing phase-shift variance with respect to sinusoidal fluctuations 
of random phases can also be defined (see Appendix B). 

Considering now changes in average light intensity, their effect on phase shift variance is easily 
derived from the property that V(0*,e) = 7 is constant, yielding 11 = [7/V(0*)]^ (see appendix B). 
Entrainment stability is affected in a more subtle manner. We find that S can be expressed in terms 
of the d-IPRC function W{(f)) and its derivatives (see Appendix B): 

r wi<p*)w"{<p*) 
[ w'{(f)*y 

where 0* satisfies eW{(f)*) = 7. If the period mismatch 7 or the second derivative W"{(f)) are zero, 
E = 1. Lower (resp., larger) values of E associated with an extended (resp., reduced) stability domain 
of the entrainment state imply 'jW"{(l)*) > (resp., < 0). Thus, robustness requires a convex PRC 
when FRP is larger than the forcing period, and a concave one in the opposite case. Again, Figs. [2] 
B and D show two d-IPRCs computed from the circadian model introduced in next section. The fact 
that their curvatures are opposite near CTO suggests that they have different robustness properties. 



(8) 



(9) 
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as we will see later. 

In this section, we have derived and discussed two simplified expressions for 11 and S which 
are valid for fiuctuations in daylight temporal profile and in average daylight intensities, respec- 
tively. This is because changes in the daylight temporal profile, unlike those in the average daylight 
intensities, tend to be uncorrelated over successive days, thus perturbing the phase-shift without 
destabilizing the entrainment state. For fiuctuations in average daylight intensities, which may per- 
sist over several days or even be permanent (corresponding to differences between two individuals), 
destabilization is the most disrupting effect. In following sections, computation of 11 (resp., S) will 
thus be carried using Eq. [8] (resp., Eq. [9]) to characterize robustness to fluctuations in daylight 
temporal proflle (resp., intensity). 



Robustness to daylight fiuctuations in a minimal model of circadian oscillator 

The criteria derived above to ensure robust entrainment in the face of daylight fluctuations hold 
for any nonlinear oscillator subjected to weak enough forcing. In this section, we check whether 
these criteria also apply to a minimal circadian oscillator model where coupling to light is generally 
non- negligible, and for which the weak forcing approximation does not necessarily hold. In most 
organisms, endogenous clock oscillations appear to rely on a core negative feedback loop, through 
which a clock gen e encodes proteins that activate (vs inactivate) its own transcriptional repressor 



vs activator) (1251 ) . The presence of delays or nonlinearities along the loop favors the emergence 



of oscillations in this autoregulatory loop (l26l . l27l ). This basic clock architecture can be captured 
in low-dimensional dynamical models such as the one originally proposed by Leloup and Goldbeter 
for Neurospora clock (28): a gene sequence is transcribed into mRNA (M), which translates into a 
protein in the cytoplasm (Pc), further translocated in the nucleus (Pn) where it inactivates the gene: 



=SpM -dp^^-hPc + hPN (10) 

=hPc-hPN 

The Michaelis-Menten-like kinetics used to describe transcription and degradation dynamics is 
required for the appearance of spontaneous oscillations. We use the following model parameters, 
which give rise to 24-hour oscillations in the dark: n = 4, sm = 2.2, Kj = 1.8, (Im = 0.84, Km = 0.5, 
sp = 0.4, dp = 1.6, kp = 0.13, ki = 0.4, ^2 = 0.45, r = 1. The effect of light is to modify one or 
several parameters during a time interval td and with gain e as in Eq. |2l To investigate the influence 
of a mismatch between the forcing period T and the FRP Tq, the latter is varied by changing the 
time constant r (Tq = 24 r). 

To illustrate our analysis, we have selected two light-coupling schemes on the basis of their IPRCs 
so that they provide us with examples of robust and non robust entrainment with respect to daylight 
fluctuations. These two schemes consist in the activation or repression of transcription by light, 
corresponding respectively to an increase or a decrease of sm during daytime. The IPRCs associated 
with transcriptional repression by light indicates an insensitivity during the time interval in which 
light is applied (Fig. [2] A) and a negative sign of the third derivative near its zero (Fig. [2] B). 
According to Eqs. [S] and [HI both properties are expected to favor entrainment robustness. The 
opposite conclusions hold for the IPRCs associated with transcriptional activation by light (Figs. [2] 
C, D), which presumably indicates poor robustness properties. 

Fig. [3] conflrms that these two forcing schemes signiflcantly differ in the robustness of the en- 
trainment state to variations in light intensity, as we anticipated from their IPRCs. In case of 
transcriptional repression by light, entrainment remains stable at various light levels (left panel of 
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Fig. [3]^), and the entrained oscillations vary little. In contrast to this, period doubling occurs when 
transcription is activated by light (right panel of Fig. |3]^). 

As predicted by the theory, these different robustness properties of the two forcing schemes when 
light level changes for a given FRP can be traced back to how the PRC and its derivative change 
when light levels are increased (Fig. O-B). In the case of transcriptional repression by light, the PRC 
varies very little in an extended neighborhood of the fixed point. In particular, ^'(0*) remains well 
within the stability boundary, whereas in the case of transcriptional activation by light, the phase- 
locked state is destabilized (V(0*) < —2) beyond a critical forcing amplitude where period-doubling 
typically occurs (Fig. |3] C). Accordingly, values of S are much lower in the case of transcriptional 
repression by light at all forcing strengths (Fig. ^D). Fig. [3] shows that these robustness properties 
do not depend on the FRP. Whatever the value of the latter, the 1:1 entrainment domain is much 
wider for transcriptional repression by light. In fact, the main effect of the period mismatch is to 
force a minimum modulation amplitude for synchronization when it it increases or decreases away 
from 0. 

These two light-coupling schemes are also associated with different robustness properties of the 
clock phase with respect to changes of the daylight temporal profile (Fig. H]). The phase shift variance 
induced by small sinusoid perturbations of different phases is much smaller (by a factor from 2 to 10) 
when light represses rather than activates the self- regulated gene (Fig. H]^). Again, this difference in 
robustness can be explained by the sentitivity or insensitivity of the PRC near the stationary phase 
shift (Fig. mi?), which ultimately depends on the existence of dead zone in the relevant time interval 
(Fig. [2] j4 and C). Fig. H] C shows that these differences in phase-shift sensitivities with respect to 
sinusoidal profile fluctuations are more or less pronounced depending on the specific value of the 
FRP To, or of the characteristic period of the fluctuations (/c = 1 or 2). 

Taken together. Figs. [3] and H] shows two light-coupling schemes that are associated with a large 
(resp. small) entrainment domain inside which there is a low (resp. large) phase-shift variability. In 
fact, the two quantities tend to be correlated, as is assessed by a systematic correlation analysis (see 
Fig. SI in the Supporting Material). 



Relating robustness to IPRC: validity of the weak coupling approximation 

The extent to which IPRC properties can account for the robustness of circadian clock with respect 
to daylight fluctuations depends on the validity of the phase approximation. In this section, we 
therefore evaluate the agreement between robustness properties derived analytically in the limit of 
weak forcing and those measured numerically when the forcing is not weak, as depicted in Fig. El 

A first indication is provided by estimating the light-stimulus amplitude beyond which the pro- 
portional relationship between phase response and light amplitude (Eq. [5]) no longer holds. In this 
respect, we found that the phase approximation can remain acceptable for modulations in excess 
50%, depending on the parameter considered and on the stimulus phase (Fig. \5\A, B). Other quan- 
tities derived from PRCs, like the phase shift 0*, the stability coefficient V'{(j)*), or the phase-shift 
and coefficient stability sensitivities 11 and S, are also expected to deviate from their estimates in the 
phase approximation. Estimation of the phase shift and the stability coefficient of the entrainment 
state is quantitatively good when parameters change by 30% between night and day (Fig. |5] C, 
D). Quantitative agreement is slightly more difficult for 11 and E (Fig. [5] i?, F), but a qualitative 
agreement can still be observed and is in fact more than sufficient to reflect a large (over 20-fold) 
dispersion of values of 11 and S between the most and least robust entrainment schemes. 

The study of the range of validity of our approach has been extended to more realistic circadian 



clock models associated with various organisms hke Drosophila (j28[ ). Ostreococcus (1l8l). cyanobacteria 
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(1291 ) or mammals (1301 ) (see Fig. S2 in the Supporting Material). The good agreement observed 
between robustness measures and their estimates using IPRC confirms that the phase approximation 
remains quahtatively valid for moderate forcing strength regardless of the model properties, thus 
providing an efficient tool to discriminate forcing schemes according to their robustness properties. 

Experimental PRCs analysis reveal robustness properties to daylight fluctuations 

We have shown that PRCs, despite of their simplicity, can provide a detailed and reliable information 
on the robustness properties of circadian entrainment. Incidentally, PRCs have been measured over 



fifty years in many organisms (1311 ). thereby offering indirect evidence of whether natural clocks are 
robust or not to daylight fluctuations. 

We have selected several PRCs measured in response to light pulses of various durations for twelve 



species (j32|-|42l) (Fig. ^A). One may assume that for type-I experimental PRCs (a to j panels), the 
phase approximation holds. Therefore, their corresponding impulse IPRCs are expected to display a 
shape similar to that of experimental PRCs if the light pulse is short enough (panels a-g of Fig. ^A) 
or can otherwise be obtained by a rectangular-pulse deconvolution (panels h-j of Fig. [6] yl). In any 
case, the impulse IPRCs display pronounced dead zones (Fig. MB), albeit of different sizes, which is 
likely to confer robustness to daylight fluctuations according to our results. 

For a quantitative assessment of this conjecture, we compare the values of 11 for the i-IPRCs 
estimated from experimental data with those obtained for many different light-coupling schemes in 
various circadian models (Fig. |6] C). Values of 11 obtained from experimental PRCs are found to be 
systematically lower than in the control case of a linear i-PRC, and are consistent with those computed 
for the most robust entrainment schemes of various circadian models. This strongly supports the 
idea that the circadian clock of these organisms have evolved robustness properties with respect to 
daylight fluctuations. 

PRCs in response to daylight-like stimuli such as 12-hour rectangular pulse also provide valuable 
information about the entrainment robustness measured by the quantity S. Unfortunately, those 
PRCs have been much less measured experimentally, and can be hardly estimated from the PRCs 
available without a precise knowledge of the light-coupling proflle. Noticeable exceptions are the 



PRCs measured in two species ( 14 ll . |42| ) that are depicted in the panels (k) and (1) of Fig. |6] A. The 
apparent change of curvature from convex to concave when the circadian phase goes from negative 
to positive seems again to indicate a robust entrainment. 



DISCUSSION 

In most living organisms, a circadian clock synchronizes internal physiological processes with cyclic 
environmental changes (e.g, in light or temperature) associated with Earth rotation. A stable phase 
relationship between internal and external time is however challenged by cellular and environmental 
variability. Part of this variability is linked to the diurnal cycle, in particular fluctuations in daylight 
perceived. In this work, we have derived explicit criteria for the robustness of circadian entrainment to 
daylight variability. Combining analytical and numerical approaches, we have shown that robustness 
depends critically on the precise shape of the IPRC, which characterizes the linear response of an 
oscillator to vanishingly small light inputs. On the one hand, the existence of an extensive range of 
entrainment is favored by a convex (resp. concave) d-IPRC (IPRC for daily pattern of light) when the 
FRP is larger (resp. lower) than the day-night driving period. On the other hand, a low dispersion of 
the phase shift between the internal oscillator and the environmental driving cycle requires a phase 
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response insensitivity of the clock during daytime, the so-called dead zone, which is manifested by 
the null interval of the i-IPRC (IPRC for light impulse) during the subjective day. 

These two PRC properties confering entrainment robustness to forcing fluctuations happen inci- 
dentally to be observed in circadian PRCs measured experimentally in many organisms, in particular 
the existence of a dead zone that seems to be universal. Although clock resetting toward a stable 
phase shift obviously requires that phase advances and delays in PRCs occur in early and late day 
respectively, the presence of a dead zone has apparently emerged as a strategy to avoid phase vari- 
ability due to daylight fluctuations and to ensure entrainment stability. The different size of dead 
zone intervals observed in various species does not weaken the robustness hypothesis but rather 
suggests the existence of light-gating mechanisms that restrict the light sensitivity of the clock to a 



certain window of the day (Il8l . 1301 ). These universal PRC properties of circadian clock may reflect 
an universal strategy to minimize the impact of fluctuations in the forcing cycle on the clock phase, 
thus suggesting a convergent functional evolution of circadian clocks. 

Trade-ofFs between multiple evolutionary goals 

Circadian clock properties, including IPRC or FRP properties, are also constrained by other adaptive 
purposes besides the need for robustness to daylight fluctuations, raising the question of whether 
trade-offs are eventually required. 

One important evolutive constraint for circadian clock is the requirement to achieve fast resetting 



after jetlag (1431 ). which depends in first approximation on the derivative of the PRC at 0* (optimal 
resetting is obtained for V'{(j)*) = — 1). For circadian oscillators that display a dead zone, the weak 
forcing approximation indicates that fast resetting requires that the dead zone must be slightly 
smaller than the daylight interval or/and that the FRP must be different from 24 hours. These two 
conditions are fully compatible with criteria for robustness to daylight fluctuations. 

The variability of the FRP of the circadian rhythm is also a prevalent source of variability sus- 
ceptible to affect clock phase (35). Daan and Pittendrigh stressed the fact that a clock characterized 
by (i) a FRP equal to 24 hours and (ii) a large dead zone of the PRC displays a phase instability for 
which a small change in the FRP causes a large change in phase shift. This result has been tempered 



by a modeling study ( ]44| ) showing that entrained circadian oscillations for the Drosophila model, 
which also display a dead zone, have a wider range of stability when the FRP is close to 24 hours. 
These apparently contradictory results can be reconciled in the light of our results when reasoning 
in terms of d-IPRCs of the circadian oscillator and depending on whether one focus on the effect 
of the FRP variability on the destabilization of the 1:1 synchronization state or on the phase-shift 
dispersion. The argument of Daan and Pittendrigh relies on the assumption that the first derivative 
of the d-IPRCs is small (i.e., the dead zone of the i-IPRC is larger than the light-coupling interval) 
whereas the argument of Kurosawa and Golbeter depends primarily on the second derivative values 
of the d-IPRC. In fact, it is likely that any circadian oscillators is characterized by an optimal value 
of the FRP whose distance from 24 hours would achieve a compromise in maintaining phase-shift 
stability and precision in face of both sources of variability in daylight intensities and FRP. 

Finally, seasonal variations of daylength is another environmental variation that the clock needs 



to adapt to by tracking dawn, dusk or both. Previous studies (1301 l45l ) have shown that the presence 
of a dead zone is an efficient mechanism to achieve seasonal tracking of dawn (resp. dusk) with a 
FRP larger (resp. smaller) than 24 hours. However, this mechanism requires that the dead zone 
is larger than the light coupling interval for some short photoperiods, which may antagonize the 
requirement for a robust synchronization in presence of daylight fluctuations. 

Thus, circadian clocks characterized by a dead zone during the subjective day and a finely tuned 
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FRP, eventually supplemented with specific light-gating mechanisms ( ll8l . |3 
an efficient trade-off between multiple and unrelated evolutionary goals 



are suitable to achieve 



Toward design principles of circadian clocks 

The concept of phase response curve was introduced and exploited to investigate theoretically and 
experimentally rhythmic and synchronization behaviors in various biological systems (|20|. l46l) . When 
oscillators are weakly perturbed, their PRC reduces to an infinitesimal analog whose linearity allows 
a more comprehensive analysis of the synchronization behavior (jlll . |22| . 1471 ) . Applying this approach 
in the context of circadian clocks allowed us to characterize the robustness of circadian entrainment 
with respect to daylight fluctuations in terms of specific characteristics of infinitesimal PRCs. This 
finding opens up a broad range of issues regarding the biophysical implementation of specific phase 
response curve properties in circadian clocks, similarly to the role of ionic currents in shaping the 
phase response curve of neuronal oscillators (j48l . |49| ). So far, many features of oscillatory behavior 
have been related to the structural basis of the underlying biochemical oscillator (jo], 26, 50-52), but 
not to phase response curve characteristics. Understanding how the phase response curve is shaped 
by various features of light coupling, including gating (30) and photo-adaptation (1531 ) . or intrinsic 
oscillator properties such as feedback, saturation or delays, would undoubtely shed light on the role 
of some key features of circadian clock designs across species. 
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Appendices 



Appendix A. General expression of robustness quantities 

The first-return map in Eq. [3] approximates the dynamics of a periodically forced non-linear oscillator 
and introduces the function namely a PRC, which is an implicit function of amplitude e and 

temporal profile L of the light forcing. We consider first that forcing properties slightly vary among 
individuals with respect to some average daylight forcing eoLo{u): 



tL{u) = eo{Lo{u) + 5tL{u)) 



'111 



where L (or Lq) is normalized with I/to L{t)dt = 1. This normalization must be preserved by 
an appropriate normalization of L that depends on the type of fiuctuations considered (appendix B). 
We can now expand the phase shift 0* and the stability coefficient d^V{(j)*) up to first order in 

6e : 

r (eo, 6e) = 0*(eo, 0) + 5e[D,,r (eo, 0)] + 0{6e^) 

d^Vir, eo, Se) = d^V{<j)*, eo, 0) + 6e[Ds,d^V{(j), eo, 0)] + 0{6e^) 

We introduce the quantities 11 and S defined as the relative variances of the phase shift and the 
stability coefficient induced by daylight fiuctuation of random amplitude: 



(12) 



'{eo,5e)-<P*{eo,0))^> 
<5e2><d^V{(p*,eo,0)2> 



(13) 



The approximation of small fiuctuation amplitudes in Eq. [12] allows to rewrite those quantities as 
following: 

n = [D,,0*(eo,O)]2 
E = [Ds,d^V{(P*,eo,0)/d^V{(P*,eo,0)]'' 



(14) 



According to the fixed point condition of the map (Eq. |3]), V{4>*{6e), 6e) = V{(f)*{0), 0) = 7. Expand- 
ing V{(j)*, 6e) up to first order in 6e leads to: 



ds.VicP*, eo, 0) + V(0*, eo, 0)dsM^o, 0) = 
n can be therefore expressed uniquely as a function of the PRC, V{(j)): 

[d^V{(l)*,eo,0)_ 

Combining Eqs. [H] and [15] allows to develop the quantity S as follows: 

d^ds,V{(j)*, eo, 0) ds,V{(j)\eo, 0)^2^(0*, eo, 0) 



(15) 



(16) 



9^r(0*,eo,O) 



(17) 
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Appendix B. Approximation of robustness quantities in the weak forcing 
limit 

Assuming that the forced periodic orbit remains in a close neighbourhood of the free- running periodic 
orbit, we can use the phase approximation, according to which: 



where 



and 



L{u)Z{u + (j))du 



L{u)Z{u + (l))du 



(19) 



(20) 



Substituting Eq. [TS]in Eqs. [TBI and [T71 allows one to express the measures of robustness introduced 
above as functions of W and W: 



n 



and 



W{<p*) 
W'{(t)*\ 

W{(t)*)W"{ 

w'{(t)*y 



(21) 
(22) 



a. Fluctuations in the average daylight intensities 

The simplest example of fluctuations in the light driving cycle corresponds to the case where the 
overall gain of the light signal fluctuates. The normalization requirement for L{u) imposes that 
L{u) = L{u) and W{(f)) = W{(f)), which allows to rewrite the sensitivity quantities of Eqs. l2T]and 
as: 



n 



7 



and 



eW'{(f)*] 

rll I A* 



1 - 



eW 



(23) 
(24) 



b. Fluctuations in daylight profiles 

Alternatively, one can also consider the case where only the temporal profile change with some 
normalized variance whereas the light intensity averaged over the day remains unchanged: 



1/td I LiuYdu = 1 
'o 



and 



TD 



Liu) du = 



(25) 



(26) 
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Note that this normahzation of L preserves the normahzation of L{u). Substituting Eqs. [T9] and [20] 
in Eq. [21] allows us to rewrite 11 as a function of Z instead of W: 



n 



Z{u + ({)*) L{u) du 
Z'iu + (t)*)L{u) du 



(27) 



One can also estimate the phase-shift sensitivity in response to sinusoidal daylight fluctuations of 
period T£,/k and phase ipk- Decomposing L{u) as a Fourier series: 



L{u,ip) = ^ lkCOs{uk/TD + 1pk) 



{2i 



and substituting Eq. [28] in Eq. [27] leads to: 

Yjk «fc cos(V'fc) + hk sin{ipk) 



Z'{u + (P*)L{u) du 



(29) 



where and bk as the kth cosine and sine Fourier coefficients of the i-IPRC truncated on the 
subinterval in which daylight perturbs the clock (usually daytime). Summing over ipk between and 
27r gives the averaged phase shift associated with a arbitrary fluctuations of zero mean {61q = 0) and 
unitary norm {J2k^k — 



n =< uiij) >^= i/2n [ u{^)d^ = V il n 

JO u 



(30) 



where 11^ are the phase-shift variances associated with sinusoidal fluctuations of period To/k : 

4 +bl 



Hk 



2W' 



k*\2 



(31) 



These quantities are for instance easily computed in the case where Z is a decreasing linear function 
on the interval of coupling, which leads to 11^ = r|,/8 k"^ tt^ for L{u) = 1 during daytime. 



Appendix C. Robustness analysis of experimental PRCs 

In this section, we briefly describe the procedure used to analyse the experimental PRCs. First, the 
fitting procedure adjusts the discrete experimental data tj, yj with j = 1, by a continuous function 
f{t). It is based on the minimization of both the fitting error and the second derivative of /: 

^1 = I E (y^ - /(^^O)' + Y Tt^''^^)]' ^32) 

j=l,N "^^ 

The ratio ^2/^1 is adjusted typically between 5 and 15 according to the data. 

For experimental PRCs that are measured using relatively short light pulse of less than one hour, 
the estimated i-IPRC, z, is assumed to be roughly equal to /. Otherwise we perform a deconvolution 
operation to extract the estimated i-IPRCs, using a genetic algorithm to find the function, z, that 
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minimize the error quantity: 



So 



fit) 



t+To 



n 2 



z{u)du 



dt 



(33) 



To estimate Il-values associated with experimental PRCs, we use Eq. 
z. 



with the estimated i-IPRCs 
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Figure Legends 

Figure [1]. 

Measures of clock robustness to daylight fluctuations. (^4) Example of two distinct daylight 
profiles {top panel) leading to different phase shifts for the entrained oscillators [bottom panel). The 
shaded area corresponds to night. [B) In the top panel, plots of the PRCs V{(j)) measuring the 
phase change induced by the daylight perturbations shown in A applied at different phases. The 
stable phase shifts 0* in the entrainment regime are solutions of V{(f)) = 7 defining . Infinitesimal 
values of A0* with respect to small forcing discrepancies defines the sensitivity measure 11. In the 
bottom panel, plot of the derivative of the two PRCs, from which the stability coefficients V'{(j)*) of 
the entrainment states and their difference can be inferred. The shaded area corresponds to stable 
entrainment, associated with —2 < V'{(j)*) < 0. Infinitesimal values of AV'{(j)*) with respect to small 
forcing discrepancies defines the sensitivity measure S. 

Figure [2]. 

Two IPRCs associated with different robustness properties and their characteristics. 

Plots of IPRCs in a circadian model where light induces either a repression {A,B) or an activation 
{C,D) of clock gene transcription. In these examples, FRP is 24 hours (7 = 0). {A,C) Plot of i- 
IPRCs ^(0) measured in response to light impulses. Whether or not they display a dead zone during 
the light-coupling interval {full line) determines the phase-shift sensitivity 11 to daylight fluctuations 
according to the Eq. [HI {B,D) Plots of d-IPRCs W{(l)) measured in response to 12-hour rectangular 
light pulses. Slight changes of 7 modify the entrainment phase shift 0* (that would correspond to 
a displacement along the full line of W). Whether or not the W"{(f)*) > (resp., < 0) {inset panel) 
when W{(j)*) > (resp., < 0) determines the sensitivity S of the entrainment stability to daylight 
fluctuations according to the Eq. O 

Figure [H 

Range of stable 1:1 entrainment in the face of light-amplitude fluctuations (^4) Time 
course over one day of the M clock component {full lines) for different light forcing amplitudes (green: 
e = 0.3; red: e = 0.6; blue: e = 0.9) which modulate negatively {left panel) or positively {right panel) 
sm {dashed line). The shaded area corresponds to night. The FRP is set to Tq = 25 hours. {B) 
Plots of V{(l)) and V'{(j)) corresponding to the different forcing strengths shown in (^4) {same color 
code). The shaded area corresponds to a stable 1:1 entrainment state. {C,D) Plots of V'{(t)*) and E 
as a function of e for Tq = 25 hours. {E) Phase diagram showing the different dynamical regimes 
that can be observed as a function of the forcing amplitude e and period mismatch 7, extrapolated 
from the existence of a stable phase shift 0* and the value of V'{(f)*) (1:1 Sync: 1:1 entrainment, 
—2 < V'{(f)*) < ; QP : quasi-periodicity, V{4>*) = 7 has no solution ; CS : complex synchronization 
and chaos, V'{(f)*) < -2). 

Figure 31 

Sensitivity of phase-shift to light-proflle fluctuations. Left and right panels correspond to 
negative and positive modulation of sm by light during daytime with e = 0.3. (^4) Time courses of 
the light-dependent parameter and of the entrained circadian oscillations over one day in presence 
of sinusoidal daylight fluctuations. The FRP To is set to 24 hours. {B) Changes in PRCs associated 
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wit fluctuations shown in A. The thick red hues indicate the range of phase shifts (7 = 0). (C) 
Corresponding values of U]/'^ as a function of 7 {full line: k = 1; dashed line: k = 2). 

Figure [5l 

Validity of the weak forcing approximation. (^4, B) Plot of normalized PRCs V{(f))/e for 
various forcing strengths ranging from e = 0.05 {black line) to e = 0.5 {red line). In the presence 
of light, sm decreases in A and increases in B. {C-F) Plots of various measures of synchronization 
and their estimates using IPRCs computed for all possible one-parameter modulations, with a light- 
modulation strength of e = 0.3 and FRP of Tq = 23.7 and 24.3 hours. {C) Phase shifts (^*. {D) 
Stability coefficient V'{(t)*). {E,F) Square roots of S and 11. Dashed lines correspond to the control 
case of an IPRC linearly decreasing during the daylight interval. 

Figure [61. 

Experimental PRCs reveal robustness properties. {A) Experimental PRC data {triangles) 



collected from many organisms (13 ll ) and fitted {full line and shaded area) according to the procedure 



described in the Appendix C. {a) Drosophila mel anog aster, 40 min pulse (1331 ): {b) Rattus albicus, 30 



min pulse (1361 ) : (c) Mus musculus, 15 min pulse (1351 ) : {d) Mesocricetus Auratus, 15 min pulse (1351 ): 



'e) Neuropora crassa, 5 min pulse (1321 ): {f) Arabidopsis Thahana, 1 hour pulse ('34); {g) Procambarus 
bouvieri, 15 min pulse (jj?! ): {h) Passer domesticus, 6 hour pulse (40); {i) Gonyaulax polyedra, 3 hour 



pulse (1391 ): {j) Paramecium Bursaria, 4 hour pulse (l38l ) {k) Nauphoeta cinerea, 12 hours pulse (1421 ): 
{I) Sarcophaga argyrostoma, 12 hour pulse (41). {B) Estimated i-IPRCs associated from the fitted 
data of the panel (a) to {j) of (^4). {C) Plots of the ranges of n^{^^-values {a-j, shaded area) computed 
from estimated i-IPRCs shown in {B) related with data shown in (^4) associated with some range of 

1 /2 

FRP around 24 hours. For comparison, plots of the ranges of H^^^-values {k, shaded area) computed 
for all coupling schemes of the computational circadian clock model (see Fig. [5]F). The horizontal 
solid line indicates the value obtained for a linear IPRC (i.e., with zero curvature). 
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